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We perform numerical simulations to investigate tidal evolution of two single-planet systems, that is, WASP-50 and GJ 1214 and 
a two-planet system CoRoT-7. The results of orbital evolution show that tidal decay and circularization may play a significant role 
in shaping their final orbits, which is related to the initial orbital data in the simulations. For GJ 1214 system, different cases of 
initial eccentricity are also considered as only an upper limit of its eccentricity (0.27) is shown, and the outcome suggests a possible 
maximum initial eccentricity (0.4) in the adopted dynamical model. Moreover, additional runs with alternative values of dissipation 
factor Q\ are earned out to explore tidal evolution for GJ 1214b, and these results further indicate that the real Q\ of GJ 1214b 
may be much larger than its typical value, which may reasonably suggest that GJ 1214b bears a present-day larger eccentricity, 
undergoing tidal circularization at a slow rate. For the CoRoT-7 system, tidal forces make two planets migrating towards their host 
star as well as producing tidal circularization, and in this process tidal effects and mutual gravitational interactions are coupled 
with each other. Various scenarios of the initial eccentricity of the outer planet have also been done to investigate final planetary 
configuration. Tidal decay arising from stellar tides may still work for each system as the eccentricity decreases to zero, and this is 
in association with the remaining lifetime of each planet used to predict its future. 
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1 Introduction 

There are 758 extrasolar planets detected up to 6 February 
2012, a large portion of which are Jupiter-like planets, as well 
as a population of super-earths with masses ranging from 1 to 
10 Earth's mass (M e ) [1]. Figure 1 shows the distribution of 
semi-major axes a and eccentricities e of all observed extra- 
solar planets by Doppler radial velocity. The eccentricities of 
the planets tend to be close to zero for very close-in orbits (a 
< 0.2), however, the eccentricities of planets are distributed 
uniformly from to 0.9 when a > 0.2, which may be ex- 
plained as a result of tidal decay and circularization in the 
exo-planetary systems. As well known, tidal dissipation is 
quite common in the planetary systems. For example, in our 
solar system, tidal effects existing in the Moon-Earth system 
may make the Moon gradually move away from the Earth, 
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and at the same time the rotation of Earth becomes a bit 
slower. In addition, tidal heating of Io (one of Galilean moons 
of Jupiter) changes the satellite's interior. Moreover, some 
planets and most satellites have been in synchronous rotation 
due to tidal dissipation [2]. 

In exo-planetary systems, the distribution of semi-major 
axes has a cut-off (0.01 AU) beyond which no planet is un- 
veiled, because very close-in planets have been destroyed by 
tidal effects. The distribution of the planetary ages and the 
semi-major axes further shows that less younger planets have 
a much smaller semi-major axes because most of them have 
been engulfed by the host star due to the stronger tides [3]. 

Herein Jupiter- like planets (super-earths) with semi-major 
axes a < 0. 1 AU, are called hot Jupiter-like planets (hot super- 
earths) [4]. However, it is difficult to elucidate how these 
planets have formed in such close-in orbits as observed ac- 
cording to current planetary formation models [5-9]. How- 
ever, a possible dynamical evolution picture may be described 
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as follows: originally, the planets form in the distant proto- 
planetary gaseous disk (several AU), then move into closer 
orbits due to the migration mechanisms such as disk migra- 
tion [10,11], the secular stellar perturbation or Kozai mech- 
anism [12] and planetary scattering (planets will be usually 
excited to high eccentricities) [13]. Once they move to a 
much closer orbit especially within 0.1 AU, tidal effects aris- 
ing from its host star then plays a vital role in the secular evo- 
lution of the planets, to shape their final orbits [6]. Finally, it 
may come into being the observed planetary configuration as 
tidal decay and circularization [14]. 

For multiple-planet systems, consisting of at least two 
planets, it is believed that gravitational interaction between 
planets is coupled with tidal process, so that this may cause 
the orbital evolution to be more complex than that of single- 
planet systems. The previous works are not well studied for 
such systems [15,16]. Therefore, we not only investigate 
tidal evolution of two single-planet systems (WASP-50 and 
GJ 1214), but also one two-planet system (CoRoT-7). When 
a planet has a nonzero eccentricity, planetary tide (arising 
from the central stellar acting on a planet) will decrease its 
semi-major axis and eccentricity until the eccentricity drops 
to zero, including current observed orbits in this process. The 
stellar tide (arising from the planet acting on the central star) 
continues to cause the planet's orbit to decay even after the 
eccentricity decreases to zero, which is associated with the 
remaining lifetime of each planet to predict its future. 
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Figure 1 The distribution of the eccentricities e and semi-major axes a [1] 
for extrasolar planets. 

2 Tidal Theory 

According to the classical tidal theory, the mean variations of 
the semi-major axis and eccentricity in a single-planet system 
are given by [16,17] 

da 4 A to 
< — >= — na^s[{\ + 23e 2 ) + le 2 D] (1) 
dt 3 



and 



de 2 c 
< — >= — nea~ 5 s[9 + ID] 
dt 3 



„ 9 k mi 5 „ 9 k\ m 5 

S=T~ °n> P = Ru 



(3) 



represent the stellar and planetary tides respectively. In the 
above equations, m,, (i = 0, 1) are the mass and radius (in 
this paper, the subscripts and 1 mean stellar and planetary 
parameters, for two-planet system, the subscripts 1 and 2 re- 
late to the inner and outer planet, respectively), a, e and n are 
semi-major axis, eccentricity and mean orbital motion, and 
ki, Qi are Love number and tidal dissipation factor, where k; 
is related to the inner tidal-effective rigidity, the distribution 
of radial density and other parameters of a planet. (2/ is a 
characterization of the tidal strength, that is larger Q, corre- 
sponding to more weaker tidal dissipation, otherwise a more 
stronger tidal interaction. 

The above tidal theory shows that a and e are nonlin- 
early coupled and cannot be treated independently as a sim- 
ple function of time. Furthermore, when the ratio of m\ /mo 
is smaller or Qo » Q\, in these cases that stellar tide can be 
ignored. Hence, the stellar tide is not taken into account in 
the planetary tidal evolution until the planetary eccentricity 
drops to be zero in this work. This indicates that the stellar 
tide predicts the remaining lifetime of the planetary system. 

3 Numerical Simulations of Tidal Evolution 

3.1 Dynamical Model 

The reference framework used here is centered at the star and 
the motion of the planet is coplanar with respect to the ref- 
erence plane. General relativity (GR) is considered for all 
planets here due to the planetary close-in orbits. Then the 
dynamical equations for a two-planet system are given as the 
following form by adding GR on the outer planet [16]: 



G{m + m{) I r 2 -ri r 2 \ 

ri = j ri + Gm 2 \ - -? r (4) 

r\ \|r 2 -ri| 3 if) 

(m +m 1 ) (ft, +f g J 
+ (ft, +f g ) + -, 



momi 
G(m + m 2 ) 



T2 + Gill] 



ri-r 2 ri 



In - r 2 | 3 r\ t 



+ (mo + m 2 ) (f tl+ f gl ) (5) 



mom 2 



mo 



where f gl and f g2 are forces acting on the inner and outer plan- 
ets generated from GR given by [16,18] 



h^^^—^h+m-vM (6) 

Gm Q m 2w Gm 2 \ A , s 1 ,ns 



(2) 



where Vi = f i and \ 2 = r 2 , are velocities of the inner and 
outer planets respectively, and c is the speed of light. 
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In the above equations, f t| and f t , are tidal forces acting on 
the inner and outer planet excited by the host star with the 
expressions of [16,19,20] 

9Gm 2 J$ 

ft . = io [ 2r '( r ' • Vl > + r i( ri x °i + Vl >] < 8) 

2Q l nir[ v 

9GmM 

^ = ~ , lo t 2 ^ • V 2> + r 2(*2 X Q 2 + V 2 )] (9) 

where G is gravitational constant, O, is the spin velocity, 
and Q'. denotes a modified tidal dissipation factor defined as 
Q'i = 3<2i/2fej, which is related to the lag time of the de- 
formation due to tidal interaction. In general, Q. = 100 is 
adopted for Earth-like planets, e.g. CoRoT-7b, but 10 5 - 10 6 
for Jupiter-like planets, for example, WASP-50b. 

In this work, we adopt a modified MERCURY6 package to 
simulate the orbital evolution by adding GR and tidal effect 
in the original codes [23], which is a general-purpose soft- 
ware package for N-body integrations and is designed to ex- 
plore the dynamical evolution of objects moving in the grav- 
itational field of a massive central body, as well to consider 
non-gravitational (e.g., cometary jet) process. Bulirsch-Stoer 
algorithm is utilized in the numerical simulations, by nature 
being slow but more accurate in the runs. The time step 
adopted here is 1/50-1/100 times the period of the innermost 
planet, with an accuracy of 10~ 12 -10~ 16 . 

3.2 WASP-50 System 

WASP-50b is a hot Jupiter discovered by transit photome- 
try in 2011 residing at a close-in orbit (0.0295 AU) around 
a G9 dwarf ( v = 11.6, 0.892 M , 0.843 R ) (Table 1) [24]. 
The current stellar activity and rotational period indicate that 
the age of the host star is about 0.8 Gyr, which is discrepant 
with the estimated age about 7 Gyr from stellar evolution. It 
may be relevant to tidal dissipation [24]. WASP-50b might 
be in the final process of undergoing tidal decay and circu- 
larization according to its current extremely tiny eccentricity 
(0.009). Considering the angular momentum of the system 
invariable throughout tidal evolution, and the initial region 
before migration due to tidal decay being more distant than 
the current location (the same situation for other systems), 
then the initial orbital elements in the numerical simulations 
are assumed as ci\ = 0.0305 AU, e\ ; n j = 0.18. Herein a 
typical tidal dissipation factor Q' { is taken as 10 s . 

Table 1 Orbital and physical parameters of WASP-50 system [24] 



Body Mass Radius Semi-major axis Eccentricity 

WASP-50 0.892 M Q 0.843 R a 

WASP-50b 1.468 Mj 1.153 Rj 0.0295 AU <0.009 



Figure 2 shows that the variations of the semi-major axis 
and eccentricity for WASP-50b. As expected, the semi-major 



axis and eccentricity decrease over secular evolution, and fi- 
nally reach present location at 0.0295 AU in a nearly-circular 
orbit within ~ 50 Myr. 
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Figure 2 Variations of the semi-major axis and eccentricity of WASP-50b 
in the tidal evolution. 

A short-period eccentric planet may possess a circular or- 
bit if its circularization timescale is shorter than the stellar age 
due to tidal evolution. The circularization timescale is a vital 
factor to reproduce the process of tidal dissipation, which can 
be described by (in a single-planet when simply considering 
the planetary tide and assuming a as constant) [25] 

e 4 / a? \i/2mi ia\ n5 
Tote = -t = -QtbrM — 7T (10) 

For Q\ = 10 5 (where gi - 26666.7 and k\ - 0.4), the circu- 
larization timescale of WASP-50 system is ~ 6.1 Myr, which 
is consistent with the numerical simulations of 9.4 Myr. 

The semi-major axis continues to suffer from tidal decay 
even its eccentricity is reduced to be about zero owing to the 
stellar tide (Figure 3). Subsequently, WASP-50b may disin- 
tegrate by the tides near the Roche limit at ~ 0.01 AU in 0.4, 
4 and 40 Gyr when taking the values of Q' as 10 6 , 10 7 and 
10 s , respectively. According to the theory, the tidal inspiral 
time of a planet is written as [26] 

°~48mWmi 
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which represents the remaining lifetime of the planet. For 
example, the very time from the current location to its Roche 
limit. As known, Q' Q plays a very important part in tidal evo- 
lution, however, Q' Q is still unknown. Hence, values of Q' Q in 
the investigated exo-planetary systems are always referred to 
those of solar system. For Jupiter-like system as WASP-50, 
when Q' Q = 10 6 , 10 7 , 10 8 , t„ = 0.29, 2.9, 29 Gyr, which is in 
good agreement with those numerical results. 



0.030 - 



Table 2 Orbital and physical parameters of GJ 1214 system [27] 




0.010 - 



Figure 3 The stellar tidal evolution of WASP-50b from the current loca- 
tion with various Q' (The dash line means the Roche limit at ~ 0.01 AU). 

3.3 GJ 1214 System 

GJ 1214b is a super-earth detected by transit method in 2009 
with a mass and radius of 6.55 M @ and 2.678 R m respectively, 
indicating that its size is between the Earth and ice giant plan- 
ets of the solar system. As for planetary structure, GJ 1214b 
resembles a planet with a composition of water surrounded 
by a hydrogen-helium envelope with only 0.05% of the total 
planetary mass [27]. Present-day observations show that GJ 
1214b is now orbiting its host star (0.153 M Q , 0.210 R Q ) at 
0.014 AU (the smallest semi-major axis in the exo-planetary 
system thus far, see Table 2) with an upper limit of eccentric- 
ity (0.27) [27]. Carter et al. (2011) [28] further constrained 
the eccentricity of GJ 1214b as 0.138 on the bases of follow- 
up light curves. However, the planetary eccentricity should 
have been deduced to be zero in a very close-in orbit due to 
tidal interaction, however this inferring is incompatible with a 
moderate eccentricity of GJ 1214b as mentioned. One possi- 
ble explanation is that there are so many unconfirmed factors 
during the process of tidal dissipation that current observa- 
tions cannot provide a reasonable constraint on GJ 1214b's 
eccentricity. In this sense, we then extensively investigate 
tidal evolution for GJ 1214b by adopting a variety of eccen- 
tricities, in attempt to understand a possible constraint on its 
initial eccentricity. Furthermore, we again perform additional 
simulations to explore the tidal evolution for this system, us- 
ing an alternative Q\ since the factor is not well determined, 
to examine a possible range of during the evolution. 



Body 


Mass 


Radius 


Semi-major axis 


Eccentricity 


GJ 1214 


0.153 Mq 


0.210 Rq 






GJ 1214b 


6.55 M e 


2.678 fi e 


0.014 AU 


<0.27 



First, to simulate the tidal evolution of various eccentric- 
ities of GJ 1214b, we assume the initial orbital elements as 
follows: «! ini = 0.0168 AU, e x ini = 0.3, 0.35, 0.4, 0.45. Tak- 
ing g'i = 100 (a typical value for terrestrial planets), Figure 
4 shows variations of the semi-major axis and eccentricity of 
GJ 1214b for the variety of the initial eccentricities, consid- 
ering both tide and GR raised by the host star. The outcomes 
indicate that final orbits are somewhat dependent on e\ ;„;, 
where the smaller the final semi-major axis (the shorter the 
circularization timescale), the larger e\ At the time GJ 
1214b finally located at 0.014 AU after experienced tidal de- 
cay, the corresponding eccentricities are deduced to be 0. 1 04, 
0.195, 0.267 and 0.330 (marked by red dots in Figure 4) re- 
spectively, of which the former three cases can meet the up- 
per limit of observed eccentricity (0.27). Hence, the possible 
largest initial eccentricity is 0.4 if the initial semi-major axis 
is 0.0168 AU. 
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Figure 4 Variations of the semi-major axis and eccentricity of GJ 1214b in 
tidal evolution. The red dots show the current locations and its corresponding 
possible eccentricities. 
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Figure 5 Tidal evolution of GJ 1214b with various Q'| ■ 

The circularization timescale is estimated to be ~ 8702.6 
yrs according to equation (10), which is in accord with our 
numerical result (~ 10,000 yrs) where Q[ = 100 (with Q\ = 
20, k\ = 0.3) in the third case {a\ m = 0.0168 AU, e\ in i = 0.4). 
Results of this run show that the final semi-major axis moves 
about 0.0136 AU on a circular orbit over ~ 5 x 10 s yrs. If it 
takes GJ 1214b so short time to circularize its orbit, the statis- 
tical results may imply that most exoplanets have circular or- 
bits, but in fact many eccentric planets have been discovered 
so far. Therefore, the circularization timescale of GJ 1214b 
does not seem to be so short possibly because the adopted 
Q\ = 100 is not reasonable, denoting that the real value of 
<2'[ would be likely to be much larger than 100. To make 
sure a potential range of Q\ , two additional runs (where Q\ = 
1000, 10000) are further investigated in the study (Figure 5). 
As seen, Q\ simply plays an important role in circularization 
timescale (where a longer circularization timescale is asso- 
cated with a larger Q' y ), but has little influence on the present- 
day eccentricity and final locations. Substantially speaking, it 
takes GJ 1214b about 0.14 Myr and 1.4 Myr to evolve into a 
circular orbit from one of possible current cases {a\ = 0.014 
AU, d = 0.267) when Q\ = 1000 and 10000, respectively. 
Therefore, it seems to be reasonable for GJ 1214b that Q' { 
is much larger than 100, indicating that this planet may own 
an alternate structure differing from that of a typical terres- 
trial planet of the solar system. In this sense, GJ 1214b is 



likely to be in a process of undergoing tidal circularization at 
a slow rate, and its circularization timescale is proportional to 
(2'j . Hence, from the simulations, we may infer that a major- 
ity of close-in eccentric exoplanets observed now have larger 
tidal dissipation factors than those of the planets in the so- 
lar system, leading to a fact that they may still bear apparent 
eccentricities far from zero. 
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Figure 6 The stellar tidal evolution of GJ 1214b with various Q' Q . 

Herein we take Q x = 100 to estimate the remaining life- 
time of GJ 1214b, in this case its final orbit locates at 0.01 18 
AU on a circular orbit. Choosing Q' Q = 10 , 10 7 , the calcu- 
lated timescale of remaining lifetime is about 28.8 and 288.0 
Gyr respectively, according to equation (11), and this is con- 
sistent with numerical results (Figure 6). The remaining life- 
time of GJ 1214b is longer than that of WASP-50b because 
of a relatively lower mass of GJ 1214. If the main-sequence 
lifetimes of stars such as the type of GJ 1214, with a mass of 
only 0. 15 M Q , are 100 Gyr by taking Q' Q = 10 6 , GJ 1214b wiU 
eventually be affected by stellar tide, however, it will survive 
throughout the stellar lifetime for Q'„ = 10 7 . 



3.4 CoRoT-7 System 

CoRoT-7b is the first super-earth discovered with a deter- 
mined radius and mass (8.0 M m , 1.58 R B ) about its host star at 
a close-in circular orbit (0.0172 AU), companied with a sec- 
ond Earth-like planet CoRoT-7c (13.6 M e , 2.39 R m ) at 0.046 
AU (Table 3) [29, 30]. Hence, the CoRoT-7 system, hosting 
two terrestrial planets, motivates a great many of researchers 
to explore its formation and dynamical evolution. 

Table 3 Orbital and physical data of CoRoT-7 system [29-31] 



Body 


Mass 


Radius 


Semi-major axis 


Eccentricity 


CoRoT-7 


0.93 M Q 


0.87 R Q 






CoRoT-7b 


8.0 M e 


1.58 R e 


0.017 AU 





CoRoT-7c 


13.6 M e 


2.39 fie 


0.046 AU 
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Figure 7 Variations of semi-major axes of CoRoT-7 system in tidal evolu- 
tion. 




Myr, next approaches a quasi-equilibrium state. The excita- 
tion of eccentricity is essential for the inner planet to undergo 
tidal decay and circularization as the tidal dissipation requires 
a nonzero eccentricity according to equation (l)-(3). After 
the relaxation progress, the orbit of the inner planet becomes 
smooth and continues to suffer from tidal decay and circular- 
ization. As a result that it stops near the current location at 
0.0175 AU in a circular orbit. Concurrently, the outer planet 
of CoRoT-7c still migrates towards the star at a much slower 
rate comparing to that of the inner planet owning to the far- 
ther distance from its host star. Although both of the system 
run towards the star at various migration rate, resulting in a 
much wider separation between them, a pair of circular orbits 
have finally formed at about 2 Myr. 
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Figure 8 The variations of eccentricities of CoRoT-7 system in tidal evo- 
lution. 

Considering the total angular momentum as invariable, the 
initial orbital elements are assumed as e\ ; n j = 0, e-i ; n ; = 0.1, 
fli ini = 0.0188 AU, a 2 m = 0.0466 AU. Figures 7 and 8 show 
the tidal evolution of two planets under consideration of both 
tide and GR {Q\ — Q' 2 = 100). The semi-major axes and ec- 
centricities of two planets are both decreasing due to tidal ef- 
fects. The initial eccentricity of CoRoT-7b is pumped up to a 
moderate value of 0. 1 due to the excitation of the outer eccen- 
tric companion, then it is relaxed within a short time of ~ 0.5 



Figure 9 Tidal evolution of CoRoT-7c in single-planet and two-planet 
cases. 

In addition, we carry out additional runs to investigate the 
tidal evolution of only a single outer planet to understand the 
outer planet's orbital decay and circularization. In this case, 
orbital shrinkage is much quicker but circularization is much 
slower than that of the two-planet case (Figure 9). As a con- 
sequence, the decreasing of ej_ is mainly produced by the tidal 
decay and circulation of the inner planet, as a result of cou- 
pled tidal and gravitational interactions rather the contribu- 
tions from its own tidal dissipation. However, to conserve the 
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total angular momentum in the planet-star system, the tidal 
decay of the outer planet is damped in presence of the inner 
planet. In theory, Mardling (2007) [32] studied the secular 
evolution of the outer planet in a two-planet system, express- 
ing the variation of e 2 as 



e 2 = -- 



e 2 
F(e 2 Y 



(12) 



where t cu - c = ei/ei is the circularization timescale of the 
inner planet in a single-planet system, and e\ is given by 
equation (12), A = (25 1 \6)(m\ I m 2 ){a\ I a 2 ) s 12 , and F(e2) = 
e\{\ - as 2 [ + ysj) 2 , a = (mi/m 2 )(ai/a 2 ) 1/2 . 
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Figure 10 Tidal evolution of CoRoT-7b with a variety of ei ;„;. 

Qualitatively, our results agree with those of Rodriguez et 
al. (2011) [16], but the resulting orbits we obtained differ 
from theirs because of the adoption of the initial orbital ele- 
ments, as the initial parameters have an important part in its 
final orbit. In this sense, three runs of simply changing e 2 ini 
= 0.1, 0.15, 0.2, respectively, are again performed to explore 
the situation of orbital evolution for a two-planet system. Fig- 
ures 10 and 1 1 show that both semi-major axes and eccentric- 
ities are influenced by e 2 ini, as a result, the final locations of 
the inner and the outer planets are at 0.0175, 0.0166, 0.0140 
AU and 0.04656, 0.04651, 0.04647 AU, respectively, and the 
changes in final semi-major axis of the outer planet is smaller 
as it is far away from the star. The final semi-major axes are 



smaller with larger e 2 ini, and the changes of their eccentric- 
ities can be ignored since all of them are deduced to be zero 
within ~ 2 Myr, implying that their final locations are mainly 
dominated by e 2 ; n i. Furthermore it provides concrete evi- 
dence that tidal effects and mutual gravitational interactions 
are coupled with each other during the secular evolution. 
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Figure 11 Tidal evolution of CoRoT-7c with a variety of ei i„;. 

Figure 12 shows the stellar tidal evolution from the current 
orbit for Q' Q = 10 s , 10 6 , 10 7 , respectively. For example, if 
taking Q' = 10 6 , CoRoT-7b will be demolished in less than 
1 Gyr, which is in concordance with the evaluated tidal in- 
spiral timescale r flJ = 0.4 Gyr. However, CoRoT-7c will still 
remain survival in the stellar lifetime due to its extremely far 
proximity from the host star by taking the above values of Q' Q . 




Figure 12 The stellar tidal evolution of CoRoT-7b with various Q' Q . 
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3.5 Summary and Discussion 

In this work, we numerically simulate two single-planet sys- 
tems including one hot Jupiter-like planetary system (WASP- 
50) and hot super-earth system (GJ 1214), and a two-planet 
Earth-like system (CoRoT-7). We may summarize the main 
results below: 

WASP-50b may arrive at its final orbit with a zero eccen- 
tricity in about 50 Myr by suffering from tidal decay and 
circularization. For GJ 1214 system, alternative initial ec- 
centricities in simulations are investigated to make clear the 
tidal evolution for the super-earth GJ 1214b. As a result that 
the final place is more closer to the star and the eccentric- 
ity of the current location is much larger with a larger initial 
eccentricity in tidal process, and an upper limit of initial ec- 
centricity (0.4) is given at the initial location (0.0168 AU). In 
addition, Q' x may have a direct influence on the circulariza- 
tion timescale of GJ 1214b, implying that Q\ of GJ 1214b 
could be probably much larger than that of terrestrial planets 
(100), so as to explain why GJ 1214b still bears an apparent 
eccentricity in such a close-in orbit. For two-planet system 
CoRoT-7, tidal evolution induces two planets to approxima- 
tively reach their current locations with zero eccentricities re- 
spectively, although both of them suffer from mutual influ- 
ences due to coupled affects of tidal and gravitational inter- 
action, which can be explained by the conservation of angular 
momentum of the planet-star system. When the eccentricity 
goes down to zero, the stellar tide then begins to work and 
causes the orbital decay of the planet until it enters the re- 
gion of Roche limit. Considering Q[ - 10 6 , the remaining 
lifetimes of WASP-50b, GJ 1214b and CoRoT-7b are about 
a few billion years in comparison with the stellar ages, how- 
ever, CoRoT-7c may survive steadily owing to relatively far- 
ther distance from its host star. 

In this paper, GR is not analyzed although it is numeri- 
cally considered in tidal evolution, because GR has no ef- 
fect on secular variations of semi-major axis and eccentricity 
[33]. Herein our numerical simulations are based on several 
assumptions: 

Firstly, tidal dissipation factor Q[ (or Q' ) is adopted from 
an equilibrium tidal theory, which is independent of tidal fre- 
quency and amplitude; secondly, Q' x is referred to those val- 
ues of the solar system. For instance, taking 10 5 and 100 for 
Jupiter-like and terrestrial planets, respectively; thirdly, the 
stellar tide arising from the planet exerting on the host star is 
ignored during tidal decay and circularization; 

Finally, the synchronization is assumed since the spin pe- 
riod of the planet is unknown and generally the synchroniza- 
tion timescale is much shorter than a secular tidal evolution. 
Thus, all assumptions may produce a deviation from an au- 
thentic dynamical evolution for the investigated systems. Cir- 
cularization timescales of some planetary systems may be un- 
reasonable when compared to the ages of planets due to the 
improper Q' v as well as the planetary ages are inconsistent 
with their stellar ages caused by the uncertain Q' . Therefore, 
new constraints of Q' { and Q' Q are expected from observations 



so as to understand tidal dissipation well. In addition, the fu- 
ture new models of planetary formation and migration are of 
assistance to restrain the starting setup of the reasonable ini- 
tial configurations adopted in numerical simulations that the 
final orbits depend on. Hence, until all above issues are well 
resolved, one complete scenario from planet formation, mi- 
gration and tidal evolution to generation of the final planetary 
configuration will help us to understand how hot Jupiter-like 
and Earth-like planets are formed. 
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